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Abstract. We study partition functions for the dimer model on families of finite graphs converging 
to infinite self-similar graphs and forming approximation sequences to certain well-known fractals. The 
graphs that we consider are provided by actions of finitely generated groups by automorphisms on rooted 
trees, and thus their edges are naturally labeled by the generators of the group. It is thus natural to 
consider weight functions on these graphs taking different values according to the labeling. We study in 
detail the well-known example of the Hanoi Towers group H^^\ closely related to the Sierpinski gasket. 
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1. Introduction 

The dimer model is widely studied in different areas of mathematics ranging from combinatorics to 
probability theory to algebraic geometry. It originated in statistical mechanics where it was introduced 
in the purpose of investigating absorption of diatomic molecules on surfaces. In particular, one wants 
to find the number of ways in which diatomic molecules, called dimers, can cover a doubly periodic 
lattice, so that each dimer covers two adjacent lattice points and no lattice point remains uncovered. 
First exact results on the dimer model in a finite rectangle of 7? were obtained by Kasteleyn [T^l [13] and 
independently Temperley and Fisher fl7| in the 60-s. A much more recent breakthrough is the solution of 
the dimer model on arbitrary planar bipartite periodic graphs by Kenyon, Okounkov and Sheffield |15j . 
We refer to the lecture notes by Kenyon [Tl] for an introduction into the dimer model. 

Let Y = {V,E) be a finite graph with the vertex set V and the edge set E. A dimer is a graph 
consisting of two vertices connected by a non-oriented edge. A dimer covering D of y is an arrangement 
of dimers on Y such that each vertex of V is the endpoint of exactly one dimer. In other words, dimer 
coverings correspond exactly to perfect matchings in Y . Let T) denote the set of dimer coverings of Y, 
and let w : E — > M+ be a weight function defined on the edge set of Y. The physical meaning of the 
weight function can be, for example, the interaction energy between the atoms in a diatomic molecule. 
We associate with each dimer covering D G 2? its weight 

W{D) H u;(e), 

eeD 

i.e., the product of the weights of the edges belonging to D. To each weight function w on y corresponds 
the Boltzmann measure on I?, /i = ^) defined as 

The normalizing constant that ensures that this is a probability measure is one of the central objects in 
the theory, it is called the partition function: 

:= W{D). 
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If the weight function is constant equal to 1, the partition function just counts all the dimer coverings on 
Y. 

For a growing sequence of finite graphs, {Yn\n, one can ask whether the limit 

log($„(u;„)) 

n^oo |y(y„)| 

exists, where Wn is the weight function on y„, and is the associated partition function. If yes, 

it is then called the thermodynamic limit. For Wn = 1, it specializes to the entropy of the absorption of 
diatomic molecules per site. 

Let us recall the method developed by Kasteleyn [T^] to compute the partition function of the dimer 
model on finite planar graphs. It consists in, given a finite graph Y ~ {V,E), constructing an anti- 
symmetric matrix such that the absolute value of its Pfaffian is the partition function of the dimer model 
on Y . Recall that the Pfaffian Pf{M) of an anti-symmetric matrix M = (TOij)ij=i,...,jv, with TV even, is 

Pf{M) := sgn{n)mp^p^ ■ ■ ■ mp^_,p^ , 

Tr£Sym(N} 

where the sum runs over all permutations tt = | ^ | such that pi < po, Ps < Pi, ■ ■ ■, 

\Pi P2 ■■■ Pn) 

Pn-1 < Pn and pi < ps < ■ ■ ■ < pn-i- One has {PJ{M)Y = det(M). 

Given an orientation on Y and a weight function w on E, consider the oriented adjacency matrix 
A = (aij)ij=i,...jy| of (y, w) with this orientation. It is of course anti-symmetric. 

Definition 1.1. A good orientation on Y is an orientation of the edges of Y such that the number of 
clockwise oriented edges around each face of Y is odd. 

Theorem 1.2 f|12j). (1) Let Y ~ (V, E) be a planar graph with a good orientation, let w be a weight 
function on E. If A is the associated oriented adjacency matrix, then 

<^{w) = \Pf{A)\. 

(2) IfY is planar, a good orientation on Y always exists. 

In this paper we apply Kasteleyn's method to study dimers partition functions on families of finite 
graphs that form approximating sequences for some well-known fractals, and on the other hand converge 
in local convergence to interesting self-similar graphs. The graphs that we consider are Schreier graphs 
of certain finitely generated groups and thus come naturally endowed with a labeling of the edges of the 
graph by the generators of the group. It is therefore natural to think about the edges with different labels 
as being of different type, and to consider weight functions on them that take different values according 
to the type of the edge. 

We now turn to Schreier graphs of self-similar groups and recall some basic facts and definitions. Let 
T be the infinite regular rooted tree of degree g, i.e., the rooted tree where each vertex has q offsprings. 
Every vertex of the n-th. level of the tree can be regarded as an element of the set X" of words of length 
n over the alphabet X = {0, 1, . . . , q — 1}; the set X" of infinite words over X can be identified with the 
set dT of infinite geodesic rays in T emanating from the root. Now let G < Aut{T) be a group acting on 
T by automorphisms, generated by a finite symmetric set S G G. Throughout the paper we will assume 
that the action of G is transitive on each level of the tree (note that any action by automorphisms is 
level-preserving) . 

Definition 1.3. The n-th Schreier graph E„ of the action of G on T, with respect to the generating set 
S, is a (labeled, oriented) graph with !/(!]„) = X"' , and edges (u,v) between vertices u and v such that u 
is moved to v by the action of some generator s € S. The edge (u,v) is then labeled by s. 

For an infinite ray ^ G dT, the orbital Schreier graph has vertex set G-^ and the edge set determined 
by the action of generators on this orbit, as above. 

The graphs E„ are Schreier graphs of stabilizers of vertices of the n-th level of the tree and the graphs 
are Schreier graphs of stabilizers of infinite rays. It is not difficult to see that the orbital Schreier 
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graphs are infinite and that the finite Schreier graphs {I]„}^]^ form a sequence of graph coverings. Finite 
Schreier graphs converge to infinite Schreier graphs in the space of rooted (labeled) graphs with local 
convergence (rooted Gromov-Hausdorff convergence [11], Chapter 3). More precisely, for an infinite ray 
^ G X" denote by ^„ the n-th prefix of the word ^. Then the sequence of rooted graphs {(!]„, ^„)} 
converges to the infinite rooted graph (S^,^ in the space X of (rooted isomorphism classes of ) rooted 
graphs endowed with the following metric: the distance between two rooted graphs (Yi,ui) and (12,^2) 
is 

Dist{{Yi,vi),{Y2,V2)) := inf <^ -;i3y^(Di,r) is isomorphic to 5^2(^2,?') \ 

Vr + 1 J 

where By{v, r) is the ball of radius r in y centered in v. 

Definition 1.4. A finitely generated group G < Aut{T) is self-similar if, for all g (z G,x X, there 
exist h Cz G,y Cz X such that 

g{xw) = yh{w), 

for all finite words w over the alphabet X . 

Self-similarity implies that G can be embedded into the wreath product Sym{q) I G, so that any 
automorphism g £ G can be represented as 

(1) .9 = ■r(.go, ■ • ■ 

where r G Sym{q) describes the action of g on the first level of T and gi G G, i = 0, ...,q — 1 is the 
restriction of g on the full subtree of T rooted at the vertex i of the first level of T (observe that any such 
subtree is isomorphic to T). Hence, if x G X and u; is a finite word over X, wc have 

g{xw) = T{x)g^{w). 

See [16| and references therein for more information about this interesting class of groups, also known as 
automata groups. 

In many cases, self-similarity of a group action allows to formulate a number of rules that allow to 
construct inductively the sequence of Schreier graphs {S,i}n>i [H US] and thus to describe inductively 
the action of the group on the n-th level of the tree. More precisely, the action of 5 G G on the n-th 
level can be represented by a permutation matrix of size g, whose entries are matrices of size q^~^ . If g 
is as in ([T]), the nonzero entries of the matrix are at position (i, t(z)) and correspond to the action of the 
restriction gi of g on the subtree of depth n — 1 rooted at i, for each j = 0, . . . , q — 1. 

In the next sections we will systematically use this description. Our idea is to define recursively an 
oriented adjacency matrix associated with the action of the generators on the n-th level, with some 
prescribed signs. The rows and columns of this matrix are indexed by the words of length n over the 
alphabet {0,l,...,q — 1}, in their lexicographic order. The signs can be interpreted as corresponding 
to a good orientation of the graph I]„ . in the sense of Kasteleyn. This allows to compute the partition 
function and the number of dimer coverings by studying the Pfaffian of this matrix. 

In this paper we compute the partition function of the dimer model on the following examples of planar 
Schreier graphs associated with self-similar actions: 

- the first Grigorchuk's group of intermediate growth (see [5] for a detailed account and further refer- 
ences); 

- the Basilica group that can be described as the iterated monodromy group of the complex polynomial 
— 1 (see |16| for connections of self-similar groups to complex dynamics); 

- the Hanoi Towers group H^^^ whose action on the ternary tree models the famous Hanoi Towers game 
on three pegs, see [9], and whose Schreier graphs are closely related to the Sierpinski gasket. Let us 
mention that counting dimers on the Schreier graphs of the Hanoi towers group H^^^ is related to the 
computation of the partition function for the Ising model on the Sierpinski triangle, via Fisher's corre- 
spondence [7], see Subsection 4.5 in our paper |5], devoted to the Ising model on the self-similar Schreier 
graphs. 

Finally we also compute the partition function of the dimer model on the (finite approximations of) 
the Sierpinski triangle. These graphs cannot be labeled so as to become Schreier graphs of a self-similar 
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group, but they are very similar to the Schreier graphs of the group H^^\ They have a few natural 
labeling of the edges in three different types, of which we describe three, and provide computations for 
two of those. 

1.1. Plan of the paper. The rest of the paper is structured as follows. In Section[2]we study the dimer 
model on the Schreier graphs associated with the action of the Grigorchuk's group and of the Basilica 
group on the rooted binary tree. Even if the model on these graphs can be easily computed directly, we 
prefer to apply the general Kastclcyn theory: the partition function at each finite level is described, the 
thermodynamic limit and the entropy arc explicitly computed. In Section [3] the dimer model is studied 
on the Schreier graphs of the Hanoi Towers group H^^K First, we follow a combinatorial approach using 
recursion and the property of self-similarity of these graphs (see Section l3.2p . A recursive description 
of the partition function is given in Theorem 13.11 The thermodynamic limit is not explicitly computed, 
although its existence is proven in two particular cases (see Proposition [33] and Proposition 13 . 5p . Then, 
the problem is studied by using Kasteleyn method (Section 13. 3^ : the Pfaffian of the oriented adjacency 
matrix is recursively investigated via the Schur complement. The description of the partition function 
that we give in Theorem 13.71 uses iterations of a rational map. In Section UJ the dimer model is studied 
on finite approximations of the well-known Sierpihski gasket: these are self-similar graphs closely related 
to the Schreier graphs of the group H^^\ Two different weight functions on the edges of these graphs 
are considered and for both of them the partition function, the thermodynamic limit and the entropy are 
computed. In Scction[5]wc perform, for the Schreier graph of H^^'' and the Sicrpinski gasket, a statistical 
analysis of the random variable defined as the number of occurrences of a fixed label in a random dimer 
covering. 

2. The partition function of the dimer model on the Schreier graphs of the 
Grigorchuk's group and of the Basilica group 

In this section we study the dimer model on two examples of Schreier graphs: the Schreier graphs of 
the Grigorchuk's group and of the Basilica group. Even if in these cases the problem can be easily solved 
combinatorially, we prefer to apply here the Kasteleyn theory because we will follow the same strategy 
in the next sections to solve the problem on more complicated graphs. 

2.1. The Grigorchuk's group. Let us start with the Grigorchuk's group: this is the self-similar group 
acting on the rooted binary tree generated by the automorphisms: 

a~e{id,id), b = e(a,c), c~e{a,d), d = e{id,b), 

where e and e are, respectively, the trivial and the non-trivial permutations in Sym{2) (observe that all 
the generators are involutions). The following substitutional rules describe how to construct recursively 
the graph I]„+i from E„, starting from the Schreier graph of the first level Ei [TJ[H]. More precisely, the 
construction consists in replacing the labeled subgraphs of I]„ on the top of the picture by new labeled 
graphs (on the bottom). 

• S • • k • • c , , d , 

uvuvuvuv 



^ ^ ^ 4 

d ^ d 

,_aQ-^ ---^-Q^, • d , , b , , c , 

lu Ou' — c — Ow lu lu Iv lu Iv lu Iv 



starting from 
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Si b,c,dC.l ~ I'O^C'' 



In the study of the dimer model on these graphs, we consider the graphs without loops. We keep the 
notation S„ for these graphs. The following pictures give examples for n = 1, 2, 3. 





In general, the Schreier graph I]„, without loops, has a linear shape and it has 2"^^ simple edges, all 
labeled by a, and 2"^^ — 1 cycles of length 2 whose edges are labeled by b, c, d. 

What we need in order to apply the Kasteleyn theory is an adjacency matrix giving a good orientation 
to E„. We start by providing the (unoriented weighted) adjacency matrix A„ of Sn, which refers to 
the graph with loops, that one can easily get by using the self-similar definition of the generators of the 
group. It is defined by putting 



and, for every n > 2, 



In-l\ I a-n-i \ / a„_i \ J I In-i 

flu 



where a„, 6„, c„, dn and /„ are matrices of size 2". If we put An = aan, Bn = bbn, C„ = cc„ and £)„ = dd„ 
then A„ is given by 

^6a„_i + ca„_i + din-i ain-i 



A,j — An + Bn + Cn + Dn — 



a/„_i bcn-i + cdn-i + db„ 



We want now to modify A„ in order to get an oriented adjacency matrix A^ for E„, corresponding to 
a good orientation in the sense of Kasteleyn. To do this, it is necessary to delete the nonzero diagonal 
entries in A„ (this is equivalent to delete loops in the graph) and to construct an anti-symmetric matrix 
whose entries coincide, up to the sign, with the corresponding entries of A„. Finally, we have to verify 
that each cycle of E„, with the orientation induced by A^, has an odd number of edges clockwise oriented. 
So let us define the matrices 



Then, for every n > 2, we put 



« '-^^ b' = \ U \ (In-. 



For each n, we put A'n ~ aa'n,Bn = bbn, C'l cc'i and Dn = dd'n- Moreover, set 
, fb + c + d \ , ^ fdln^i ^ \ f 

where, for every n > 1, the matrix J„ is obtained from J„ with the following substitutions: 

b ^ d ci-^ b d 1-^ c. 
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Define 

■ a 
-a 



A'l = A[ +B[+C[+ D[ - Ji = 
and, for eacli n > 2, 



K = K + B' + C' + D' - J„ = 



Note tliat tlie matrix J„ is introduced to erase tlie nonzero diagonal entries of A„, corresponding to loops. 

Proposition 2.1. The matrix induces a good orientation on the Schreier graph E„ of the Grigorchuk's 
group. 

Proof. It is easy to show by induction that AJ^ is anti-symmetric and that each entry of A^ coincides, 
up to the sign, with the corresponding entry of the adjacency matrix A„ of E„, where loops are deleted. 
Finally, we know that all cycles in the Schreier graph have length 2 and this ensures that each cycle has 
a good orientation in the sense of Kasteleyn. □ 

Theorem 2.2. The partition function of the dimer model on the Schreier graph I]„ of the Grigorchuk's 
group is 

<i>„(a, 6, c, d) = 

Proof. It is easy to check, by using the self-similar definition of the generators of the group, that 
and 

a(l") = 01"-^ 5(1") ^ c(l") = d{r) = 1". 

This implies that the vertices 1"~^0 and 1" are the (only) two vertices of degree 1 of E„, for each n. This 
allows us to easily compute det(A^) by an iterated application of the Laplace expansion. We begin from 
the element a at the entry (2"~^,2"), which is the only nonzero element of the column 2". So we can 
"burn" the row 2"~^ and the column 2". Similarly, row 2" and column 2"~^ can be deleted and a second 
factor a appears in det(AJ^). With these deletions, we have "deleted"in the graph all edges going to and 
coming from the vertex 01"~^ (corresponding to the row and column 2"~^). So the vertex 001""^ (which 
is adjacent to 01"~^ in S„) has now degree 1 and on the lines corresponding to it there is just a letter 
a (or —a) corresponding to the edge joining it to 101""^. Hence, the Laplace expansion can be applied 
again with respect to this element, and so on. Observe that each simple edge labeled a of S„ contributes 
to det(A^). The assertion follows since the number of simple edges is 2"~^. □ 

Corollary 2.3. The thermodynamic limit is ^loga. In particular, the entropy of absorption of diatomic 
molecules per site is zero. 

Proof. A direct computation gives 

log($„(a,6,c,d)) log($„(a,6,c,d)) 1 
lim ; — = lim = - log a. 

n^+oo |V'(i;„)| «^+oo 2" 2 

By putting a = 1, we get the entropy. □ 

2.2. The Basilica group. The Basilica group [TD] is the self-similar group generated by the automor- 
phisms: 

a — e{b, id)., b — e(a, id). 

It acts level-transitively on the binary tree, and the following substitutional rules [S] allow to construct 
inductively E„+i from E„, 
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Ou 



Ov 



a 

Ou 



Ov 



starting with the Schrcicr graph Si on the first level. 



We consider here the dimer model on the Schreier graphs of the Basilica group without loops, as in the 
following pictures, for n = 1, 5. 





b b 
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It follows from the substitutional rules described above that each E„ is a cactus, (i.e., a separable graph 
whose blocks are either cycles or single edges), and that the maximal length of a cycle in S„ is \^] . We 
further compute the number of cycles in E„, that will be needed later. Denote by a* the number of cycles 
of length j labeled by a in and, similarly, denote by 6* the number of cycles of length j labeled by b 
in S,;. 



Proposition 2.4. For any n > 4 consider the Schreier graph I]„ of the Basilica group. For each k > 1, 
the number of cycles of length 1^ labeled by a is 



2«-2fe-i 1 < A: < I - 1 



and the number of cycles of length 2*^ labeled by b is 
2 k = ^ , forn odd, 



1 



n-1 
2 

L _ n+1 

K — 



02k 




ft — -2 



1 < fc < t 



for n even 



for n even. 



Proof. The recursive formulae for the generators imply that, for each n > 3, one has 



and 



"2 



-)ri-2 



and in general a2k = 0-2 '^^^ ^2*= ^ ^2 ^^'^ ■ particular, for each n > 4, the number of 2-cycles 

labeled by a is 2"^^^ and the number of 2-cycles labeled by 6 is 2"^^. More generally, the number of 
cycles of length 2*^ is given by 

n rin—2k—l in c\n—2k 

CL2k — ^ , 02k — ^ , 

where the last equality is true if n — 2fc + 2 > 4, i.e., for fc < § — 1- Finally, for n odd, there is only 
one cycle of length labeled by b and there are four cycles of length , two of them labeled by a 
and two labeled by 6; for n even, there are three cycles of length 2^ , two of them labeled by b and one 
labeled by a. □ 

Corollary 2.5. For each n > 4, the number of cycles labeled by a in S„ is 

( 2ri+2 n odd, 
n even 



Zl+1 

3 



and the number of cycles labeled by b in I]„ is 




n odd, 
n even. 
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The total number of cycles of length > 2 is 2" ^ + 1 and the total number of edges, without loops, is 
3 • 2"-i. 

Also in this case wc construct an adjacency matrix associated with a good orientation of S„, in 
the sense of Kasteleyn. We first present the (unorientcd weighted) adjacency matrix A„ of the Schreier 
graph of the Basihca group. Define the matrices 

ai = flj^^ = ^ , and hi = b^^ ^ (l 
Then, for every n > 2, we put 

In-lJ ' ~ In-lJ ' " " [in^l j ' ~ [a-\ 

If we put An ~ aon, A^^ = Q-Q-n — bbn and B^^ = bb^^ , then A„ is given by 

We modify now A„ in order to get the oriented adjacency matrix A^. To do this, we need to delete the 
nonzero diagonal entries and to construct an anti-symmetric matrix whose entries are equal, up to the 
sign, to the corresponding entries of A,i. Finally, we have to check that each elementary cycle of I]„, with 
the orientation induced by A^, has an odd number of edges clockwise oriented. We define the matrices 

I) ?) ' = ("o^ -l) ^^'i ^'i = = (-1 I, 

Then, for every n > 2, we put 



In^lJ' ^" V -In-lJ' " \-In-l J' " 

(Observe that here the exponent —1 is just a notation and it does not correspond to the inverse in 
algebraic sense.) Put A'^^ = aa'j^,A'~^ = aa'^^,B'^ — bb'^ and B'~^ = bb'~^. Then 

26^ 



A[ = A[ + A[-^ + B[ + B[-^ 

and, for each n > 2, 

a: ^ al + a'-^ + bl + b'-^ 



-2b 



6K-\-/„_i) 
Proposition 2.6. induces a good orientation on the Schreier graph S]„ of the Basilica group. 

Proof. It is easy to show by induction that A^ is anti-symmetric and that each entry of A^ coincides, 
up to the sign, with the corresponding entry of the adjacency matrix A„ of E„, where loops arc deleted. 

We also prove the assertion about the orientation by induction. For n = 1,2 we have A'^ = ^ 26 0^ 

/ 2a 25 0\ 

-2a 25 I , . , 

which correspond to 



and A 2 



-26 

V -26 oy 




10 •: 



Now look at the second block 6(ajj„]^ + In-i) — b 
corresponds to the 2-eycles 
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01m» 




•llu 



which come from the a-loops of E„_i centered at lu and so they are 2"~^. The block b{bj^_2 + 
corresponds to the 6-cycles of length 2*^ > 4. These cycles come from the a-cycles of level n — 1 but 
they have double length. In particular, &&^_2 corresponds to the 6-cycles at level n — 2 (well oriented by 
induction), that correspond to the a-cycles at level n — 1 with the same good orientation given by the 
substitutional rule 



• - 

u 



Ou 



Such a cycle labeled a with vertices ui,M2, • ■ • ,U2k-i (of length 2*^^^ > 2) at level n — 1 gives rise to a 
6-cycle of length 2'^ in S„ following the third substitutional rule. In this new cycle, by induction, there 
is an odd number of clockwise oriented edges of type 



• - 

Ou,: 



— • 

lu. 



i+l 



All the other edges have the same orientation (given by the matrix 6/„_2). Since these edges are in even 
number, this implies that the 6-cycle is well oriented. A similar argument can be developed for edges 
labeled by a and this completes the proof. □ 



Theorem 2.7. The partition function of the dimer model on the Schreier graph S„ of the Basilica group 
is 

n odd, 



$„(a,6) 



'2^6^" 



Proof. For small n the assertion can be directly verified. Suppose now n > 5. Observe that we have 
det(AJ^) = 6^" (det(ajj_;^ + /,i„i))^, since the matrices a',^_i + /„_i and a^~_\ — /„„! have the same 

determinant. Let us prove by induction on n that, for every n > 5, {det{a'^_i + = 2^" ^ • 2^' , 

where V is the number of cycles labeled by b in S„ having length greater than 2. One can verify by direct 
computation that det(A[;) = 2^^ det(A^) = 2^''*. Now 



dct(A:,) 



(det(a:,_i+/„_i))^ 



bl 



Vi-2 + 









(det(6:,_2 + ^-2))' 



-,2"- 



— In-3 



In- 



(det(a:,_3 + /„-3))' 



-)2;" 



where the last equality follows by induction and I" is the number of 6-cycles in S,i_2 having length greater 
than 2. Now observe that I" is also equal to the number of a-cycles of length greater than 2 in but 
also to the number of 6-cycles of length greater than 4 in I]„. We already proved that 6J = 62 = 2"~^ 



and so 2^" = 2'^''" 
in Proposition 12.41 



Similarly 2^ 



2262 xhen one gets the assertion by using computations made 



□ 



Corollary 2.8. The thermodynamic limit is i log 2 + ^ log 6. In particular, the entropy of absorption of 
diatomic molecules per site is ^log2. 
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3. The dimer model on the Schreier graphs of the Hanoi Towers group H^^'> 



We present here a more sophisticated example of dimers computation on Schreier graphs - the Schreier 
graphs of the action of the Hanoi Towers group i/^'^^ on the rooted ternary tree. 



3T. The Schreier graphs. The group H'^^'> is generated by the automorphisms of the ternary rooted 
tree having the following self-similar form [3]: 

a {01){id,id,a) b = (02) {id, b, id) c = {12){c,id,id), 

where (01), (02) and (12) are transpositions in Sym{3). Observe that a, b, c arc involutions. The associated 
Schreier graphs are self-similar in the sense of [12], that is, each ^n+i contains three copies of glued 
together by three edges. These graphs can be recursively constructed via the following substitutional 
rules [5]: 

llu 



Rule I 




OOu " lOu ^ I2u ^ 22u 



22u 



Rule II 




o'ou ^ 20m " 2'lu llu 



Ou 



OOu 



lu 



Rule HI 



Rule IV 



llu 

b 



2u 



22u 



Rule V 



Ov 



OOv 



Iv 



llv 



2v 



22v 



The starting point is the Schreier graph Ei of the first level. 
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In fact, the substitutional rules determine not only the graphs I]„ but the graphs together with a particular 
embedding in the plane. Throughout the paper we will always consider the graphs embedded in the plane, 
as drawn on the Figures, up to translations. Observe that, for each n > 1, the graph I]„ has three loops, 
at the vertices 0", 1" and 2", labeled by c, 6 and a, respectively. Moreover, these are the only loops in 
S„. Since the number of vertices of the Schreier graph I]„ is 3" (and so an odd number), we let a dimer 
covering of S„ cover either zero or two outmost vertices: we will consider covered by a loop the vertices 
not covered by any dimer. For this reason we do not erase the loops in this example. 
The two subsections below correspond to the calculation of the dimers on Hanoi Schreier graphs by two 
different methods: combinatorial fSection l3.2[) . using the self-similar structure of the graph; and via the 
Kasteleyn theory fSeetion l3.3|) . using self-similarity of the group H'^^^ in the construction of the oriented 
adjacency matrix. 




At level 2, we have eight possible dimer configurations: 





More generally, we will say that a dimer covering is of type I if it contains all the three loops, of type 
II if it only contains the leftmost loop (at vertex 0"), of type III if it only contains the upmost loop (at 
vertex 1") and of type IV if it only contains the rightmost loop (at vertex 2"). 

For S„, n > 1, let us denote by $^(a, b, c) the partition function of the dimcr coverings of type i, for 
i = I, II, III, IV, so that $„ = <I>^ + + <^l^^ + <tl^ . In what follows we omit the variables a, 6, c in 
the partition functions. 

Theorem 3.1. The Junctions {'^n\, n=\, ...,4, satisfy the system of equations 



^n+l V "7 abc ' n ^' 

, = (<^^^f ■ i + 

n+1 V "'/ c 71 n n 

\ 71 J b 71^ n ^71 

^ n+l \ 71 } a 71 n n 



with the initial conditions 



'$( = abc 



Proof. We prove the assertion by induction on n. The initial conditions can be easily verified. The 
induction step follows from the substitutional rules. More precisely, in I]„ we have a copy Tq of E„_i 
reflected with respect to the bisector of the angle with vertex 0"""'^; a copy Ti of E„_i reflected with 
respect to the bisector of the angle with vertex 1"~^; a copy T2 of S„_i reflected with respect to the 
bisector of the angle with vertex 2"~^. 




Using this information, let us analyze the dimer coverings of I]„ as constructed from dimer coverings of 
Ti, i = 0,1, 2, that can be in turn interpreted as dimer coverings of 

First suppose that the dimer covering of E„ contains only one loop, and without loss of generality 
assume it is at 0". There are then two possible cases. 

• The copy of S]„_i corresponding to To was covered using three loops (type I). Then the covering 
of E„ must cover the edges connecting Tq to Ti and T2 , (labeled a and b respectively) . So in the 
copy of corresponding to Ti we had a dimer covering with only a loop in A (type IV, by 

reflection) and in T2 a covering with only a loop in D (type III, by reflection). These coverings 
cover vertices E and F and so the edge labeled c joining E and F does not belong to the cover 
of E„. We describe this situation as: 




I III 

The copy of S„_i corresponding to Tq was covered with only one loop in 0"^^ (type II), so that 
the vertices B and C arc covered. This implies that the edges joining A, B and C, D, labeled 
a and b respectively are not covered in E„. So in the copy of E„_i corresponding to Ti there 
was no loop at (or A), which implies that in it only the loop at F was covered (type II, by 

reflection). Similarly the covering of the copy of E„_i corresponding to T2 could only contain 
a loop in E (type II, by reflection). Consequently, in E„, the edge joining E and F (labeled c) 
must belong to the covering. Schematically this situation can be described as follows. 

II 




Now suppose that the covering of I]„ contains loops at 0", 1" , 2". There arc again two possible cases. 

• We have on Tq a dimer covering with three loops (type I) , so that the dimer covering of E„ must 
use the edges joining the copy Tq to the copies Ti and T2. So in the copy of E„_i corresponding 
to Ti (and similarly for T2) we necessarily have a covering with three loops (type I). 




I I 

We have on Tq a dimer covering with only one loop in 0"~^ (type II), so that the vertices B and 
C are covered. This implies that the edges joining A,B and C,D are not covered in E„. So in 



COUNTING DIMER COVERINGS ON SELF-SIMILAR SCHREIER GRAPHS 



15 



Ti we cannot have a loop at A, which imphes that the corresponding copy of was covered 

with only a loop in 1"^^ (type III). Similarly T2 was covered with only a loop in 2"^^ (type IV). 



Ill 



II IV 

The claim follows. □ 
Corollary 3.2. For each n > 1, the number of dimer coverings of type I, II, III, IV of the Schreier graph 

3"-l-l 

Tin is the same and it is equal to 2 2 . Hence, the number of dimer coverings of the Schreier graph 
I]„ is equal to 2 2^ . The entropy of absorption of diatomic molecules per site is ^ log 2. 



Proof. By construction in the proof of Theorem 13.11 we see that the number of configurations of type I of 
Tin is given by 2h'^ , where h is (by inductive hypothesis) the common value of the number of configurations 

of type I, II, III, IV of equal to 2 2 . So the number of configurations of type I in I]„ is equal 

to 

3(3"~2_i) 3"-l_i 

2 • 2 2 = 2 — 2 — . 

Clearly the same count holds for the coverings of type II, III and IV of S„, and this completes the 
proof. □ 

Remark 3.3. Analogues of Theorem 3.1. can be deduced for the dimers partition function on Schreier 
graphs of any self-similar (automata) group with bounded activity (generated by a bounded automaton). 
Indeed, V. Nekrashevych in [TB] introduces an inductive procedure (that he calls "inflation" ) that produces 
a sequence of graphs that differ from the Schreier graphs only in a bounded number of edges (see also 
[3], where this construction is used to study growth of infinite Schreier graphs). This inductive procedure 
allows to describe the partition function for the dimer model on them by writing a system of recursive 
equations as in 

Unfortunately, we were not able to find an explicit solution of the system We will come back to 
these equations in Section [5] where we study some statistics for the dimer coverings on S„. Meanwhile, 
in the next Subsection 13. 3[ we will attempt to compute the partition function in a different way, using 
the Kasteleyn theory. 

In the rest of this subsection we present the solution of the system ^ in the particular case, when all 
the weights on the edges are put to be equal, and deduce the thermodynamic limit. 

Proposition 3.4. The partition function for a ^ b = c is 

$„(a, a, a) = 2~ 2 - . a^~^ (a + 3). 
In this case, the thermodynamic limit is ^log2 -f- ^loga. 

Proof. By putting a = b = c, the system ([2]) reduces to 
(3) 



3 

n+1 



n+l 

with initial conditions ${ = and ~ a? , since = = ■ One can prove by induction that 



a, so 



that the first equation in ([3]) becomes 



n+l „3 : 
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giving 



$„ = 2 2 . a 2 

;^TT n """^-! 3" + l 



The partition function is then obtained as 

$„ $^ + 3$,V = 2^^^ • (a + 3). 

The thermodynamic limit is 

hm = - log 2 H — log a. 

Finally, by putting a = 1, we get the entropy of absorption of diatomic molecules as in Corollary [321 D 

Next, we deduce the existence of the thermodynamic limit of the function $„(1, 1, c), for c > 1. Observe 
that $„(l,l,c) = $^(l,l,c) + $^^(l,l,c) + 2$,V^(l,l,c), since $^^^(1,1, c) ^ $^^(1,1, c). A similar 
argument holds for the functions $„(a, 1, 1) and $n(l, !)■ 

Proposition 3.5. For every c > 1, the thermodynamic limit lim„_^oo gi"'"'^''''*'* exists. 

Proof. It is clear that, for every c > 1, the sequence £„ := '"^^^gl^'"^''^^^ is positive. We claim that e„ is 
decreasing. We have: 

£«+! ^ log($„ + i(l,l,c)) 

En 3-l0g($„(l,l,c)) 

3.1og($„(l,l,c)) + log( ^^;YM:cpO 
3-log($„(l,l,c)) 

f <l>„+i(i4,c) ^ 



= 1 



3-log($„(l,l,c))- 



Since log($„(l, l,c)) > for every c > 1, it suffices to prove that "^"Y^ j'^y? is less or equal to 1 for 
every c > 1. 



<i>n+l(l,l,c) ^ + 1,C) + 2ci>,y/,(l, l,c) 

$„(1, 1, C)3 ($^(1, 1, c) + $,^(1, 1, C) + 2$^"(1, 1, C))3 



($^ + $// + 2ci>^"f 

< 1 

This implies the existence of lim„^oo '"^'''^gl"'"'"'"''^'''' . □ 

3.3. Partition function by Kasteleyn method. In order to define a matrix inducing a good orienta- 
tion, in the sense of Kasteleyn, on the Schreier graphs of H^^\ we introduce the matrices 



ai = -1 , fei = 1 , ci 
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Then, for every n even, we put 












n-1 







bn 








an-i/ 











In- 





bn-1 





rin- 


-1 











-In- 





bn-l 





/n-1 









and for every n > 1 odd, we put 

In-l 

-/n-1 
ttn-l^ 

where a„,6„,c„ and /„ are square matrices of size 3". Now we put An = aa,i,/3„ = bbn and C„ = cc„ 
for each n > 1 and define A„ = A„ + Bn + C„ , so that 



Cn-1 














— /n-1 





/n-1 





Cn-1 














/n-1 





— /n-1 






and, for each n > 1, 







a 


-b 




—a 


b 






[b 


— c 


a 





f Cn-1 


-a/n-1 


bin-l 


An= 1 


aIn-1 


Bn-1 


-c/n-1 




\-bIn-l 


c/„-i 


^n-1 







( Cn-1 


a/n-1 


-fo/n-1 


for n even. 




-a/n-1 


/3n-l 


c/„-i 






^ bIn-1 


-c/n-1 


An-1 



for n odd. 



We want to prove that, for each n > 1, the oriented adjacency matrix A„ induces a good orientation on 
S„. Then we will apply Kasteleyn theory to get the partition function of the dimer model on E„. One 
can easily verify that also in this case the entries of the matrix A„ coincide, in absolute value, with the 
entries of the (unoriented weighted) adjacency matrix of the Schreier graphs of the group. The problem 
related to loops and their orientation will be discussed later. The figures below describe the orientation 
induced on Si and S2 by the matrices Ai and A2, respectively. 

b 





Proposition 3.6. For each n > 1, the matrix A„ induces a 



orientation on 



Proof. Observe that in Si the sequence of labels a, 6, c appears in anticlockwise order. Following the 
substitutional rules, we deduce that for every n odd we can read in each elementary triangle the sequence 
a, 6, c in anticlockwise order. On the other hand, for n even, the occurrences of a, b, c in each elementary 
triangle of E„ follow a clockwise order. We prove our claim by induction on n. For n = I, the matrix Ai 
induces on Ei the orientation shown in the picture above, so that the assertion is true for n = 1. Now 
observe that, for every n odd, the blocks ±a/„_i, ±6/„_i, ±c/„_i in A„ ensure that each elementary 
triangle in E„ has the same orientation given by 
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lu 



a/ \c 

1 • 

Ou ° 2u 

For n even, the sequence a,b,c is clockwise and the blocks ±a/,i_i, ±6/„_i, ±c/„_i in A„ ensure that 
the orientation induced on the edges is clockwise as the following picture shows: 

2u 
• 

a • 

Ou lu 

So we conclude that for every n all the elementary triangles of I]„ are clockwise oriented. Now construct 
the graph S„_|_i from S„ and suppose n odd (the same proof works in the case n even). Rule I gives 




OOu ° 10m ^ 12u ^ 22u 



In order to understand why the edges not belonging to an elementary triangle have this orientation, we 
observe that the edge 

• ^ • 

20u 21m 

has the same orientation as the edge 

• a • 

Om 1m 

since the entry (20m, 21u) of the matrix A„+i is the same as the entry (Om, 1m) of the matrix A„. Similarly 
for the other two edges joining vertices 01m, 02u and 10m, 12m. This implies that each elementary hexagon 
has a good orientation. Now note that in T,n+i we have 3"~* elementary cycles of length 3-2*, for each 
i = 0, 1, . . . ,n. We already know that in all the elementary triangles and hexagons have a good 

orientation. Observe that each cycle in I]„ having length fc = 3 • 2™, with m > 1, gives rise in S„+i to a 
cycle of length 2k. In this new cycle of En+i, k edges join vertices starting with the same letter and keep 
the same orientation as in S„ (so they are well oriented by induction); the remaining k edges belong to 
elementary triangles and have the form 



xu xu 
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where x and x,x & {0, 1, 2}, u e I]„. Since the last k edges belong to elementary triangles, they are 
oriented in the same direction and, since k is even, they give a good orientation to the cycle. The same 
argument works for each elementary cycle and so the proof is completed. □ 



The matrix A„ cannot be directly used to find the partition function because it is not anti-symmetric 
(there are three nonzero entries in the diagonal corresponding to loops) and it is of odd size. Let r„.c be 
the matrix obtained from A„ by deleting the row and the column indexed by 0" and where the entries 
(1", 1") and (2", 2") are replaced by 0, so that the partition function of dimer coverings of type II is given 
by c-\/det(r„_c). Similarly, we define r„_f,,r„_a for dimer coverings of type III, IV, respectively. Now let 
A„ be the matrix obtained from A„ by deleting the three rows and the three columns indexed by 0", 1" 
and 2", so that the partition function of the dimer coverings of type I is abc^/ del A„. This gives 

(4) $„(a, 6, c) = cy^det(r„,c) + 5y^dct(r„,b) + a^dci{Tn,a) + abc^ dct A„. 

In order to compute det(r„.c) (the case of r„^b,r„.a and A„ is analogous), we put 












6'i 



Then for every n > 1 odd wc put 






ain-i 






aIn-1 















a'n-ij 





and for every n even we put 

























^-ra-l 



b' = 



Finally, set 



Then for every n > 1 odd wc put 








"-'n-1 














-&/n- 





b'n-1 





^n-l 














bln- 





b'n-l 





bln- 


1 








^n-1 





-blU 








-c/n-1 





c/n-1 








c/n-1 





-c/n-1 




Cn 



and for every n even we put 



<-l 



"-'n-1 






with 



/I 




Finally, let A„ = An + Bn + C„ for each n> 1, so that 






-hi" 

''-'n-1 






0/ 




Cn-1 


^I'n-l 


-bll-1 


-<-l 


b'n-1 


c/„-i 


biLi 


-c/„-i 


a'n-l 



Cn 







/ c;,_i 


-<-i 


bl'n-1 


for n odd. 


A„- 


<-i 


b'n-l 


-c/n- 






\-blLi 


c/„-i 


a'n-1 



for n even. 
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The introduction of the matrices guarantees that det(A„) = det(r„_c), since we have performed all the 
necessary cancellations in A„. Geometrically this corresponds to erasing the loops rooted at the vertices 
1" and 2" and the edges connecting the vertex 0" to the rest of I]„. 
Next, we define a rational function F : — as follows 

/ Xixl + X2X3X5Xe X2xl + XiXsXiXe XsX^ + XiX2X4X5\ 
F{xi,X2,X3,X4:,X5,Xq) Xi,X2,X3, ■ , ■ , ■ . 

\ X1X2X3 + X4X5XQ X1X2X3 + X4X5Xe X1X2X3 + x4X5Xe J 

Denote F^'^^ (x) the fc-th iteration of the function F, and Fi the i-th projection of F so that 

Fix) = iF,ix),...,Fe{x)). 

Set 

F('=)(a,6,c,a,6,c) = (a, 6, c, a^, c^'^'). 

Theorem 3.7. For each n > 3, the partition function <f>„(a, &, c) of the dimer model on the Schreier 
graph I]„ of the Hanoi Tower group IS 

71 — 3 



$„(a, b,c) = H (abc + aC^^^C^^cf^)) (a6c (a("-2)j,(«-2) + a("-2)^(„-2) ^ ^(n-2)^(„-2) 

fc=0 

+ a("-2)^(„-2)^(n-2) ^ ^^^j ^ a2(^(„-2))3 ^ ;,2(^(„-2))3 ^ ^2(^(„-2))3^ 



Proof. Wc explicitly analyze the case of dimer coverings of type II. It follows from the discussion above 
that $f/(a, 6,c) = cWdet(A„). More precisely, the factor c corresponds to the label of the loop at 0", 



and the factor ydet(A„) is the absolute value of the Pfaffian of the oriented adjacency matrix of the 

graph obtained from I]„ by deleting the edges connecting the vertex 0" to the rest of the graph and the 
loops rooted at 1" and 2". The cases corresponding to coverings of type III and IV are analogous. If we 
expand twice the matrix A„ using the recursion formula and perform the permutations (17) and (58) for 
both rows and columns, we get the matrix (for n odd) 



(Mil Mi2\ 
VAf21 M22J 












~cI„-2 


— aIn-2 





bI°-2 














-cIn-2 





—bI„-2 








aln~2 








cIn-2 











aI„-2 








-bln^2 


cIn-2 














bin -2 


-<-2 








aIn-2 


bIn-2 

















-c/„-2 













-bin-. 2 














cIn-2 


-ft/-" 
''-'n-2 



























—aI„-2 








cI„-2 








b'n-2 











bl„~2 








-c/„-2 








<-2 , 



V 

Note that each entry is a square matrix of size 3"~^. Hence, the Schur complement formula gives 



(5) 



det(A„) 



det(Afii) • det(M22 - MaiAff^^A'/is) 

c: 



(2a6c) 



2-3" 



2abc -'"-2 



n-2 

a'*+b^c^ TO y 
2abc "n-2 

b^+a^c'^ tO c-^+a^b^ 
2abc -'ra-2 2abc 



b'^+a^c'^ tO 
2ab^ Ji-2 



c^+a^b-' 
2abc 



^n-2 



The matrix obtained in ([5]) has the same shape as A„_i, so we can use recursion by defining 

-F4(xi, 



Akixi, 



F4{xi, . . . ,a;6)/°_i 
^-i^5(a;i, . . . ,a;6)/°_i 



Fe{xi, ...,xg 



-Fg{xi,. ..,X6) 



''k-l 
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Hence, ^ becomes 

det(A„) = {2abcf'^"''' ■ det(An-i{a,b,c,a,b,c)) 

= (2a5c)2-3""' (abc + a'^^^b^^^ c'^^^)^'^"'' ■ dot (a„_2 (F'^^^ (a, b, c, a, b, c))) 

= ]J (a6c + a(''^)6('^'c('^-)) ' • dct (AsIF^"-^) (a, 6, c, a, &, c))) 

/c=0 

A similar recurrence holds for coverings of type III and IV. For dimcr coverings of type I, by using the 
Schur complement again, we get 

det(A„) = n (a6c + a«6«c('=))'''" ' ' • (a6c + a("-2)^(„-2)^(„-2)^ ' ^ 

Then we use ^ and the proof is completed. □ 

Remark 3.8. The proof above, with a = 6 = c = 1, gives the number of dimers coverings of I]„ that we 
had already computed to be 2 2^ in Corollary 13.21 

4. The dimer model on the Sierpinski gasket 

In this section we study the dimer model on a sequence of graphs {r„}„>i forming finite approximations 
to the well-known Sierpinski gasket. The local limit of these graphs is an infinite graph known as the 
infinite Sierpinski triangle. The graphs r„ are not Schreier graphs of a self-similar group. However, they 
are self-similar in the sense of [12], and their structure is very similar to that of the Schreier graphs 
S„ of the group H^^\ studied in the previous section. More precisely, one can obtain r„ from S„ by 
contracting the edges joining two different elementary triangles. 

One can think of a few natural ways to label the edges of these graphs with weights of three types. The 
one that springs first into mind is the "directional" weight, where the edges are labeled a, b, c according to 
their direction in the graph drawn on the plane, see the picture of T3 with the directional labeling below. 
Note that the labeled graph r„ is obtained from the labeled graph r„_i by taking three translated copies 
of the latter (and identifying three pairs of corners, see the picture). 




The dimer model on r„ with this labeling was previously studied in |18j : the authors wrote down a 
recursion between levels n and n + 1, obtaining a system of equations involving the partition functions, 
but did not arrive at an explicit solution. Unfortunately, we were not able to compute the generating 
function of the dimer covers corresponding to this "directional" weight function either. Below we describe 
two other natural labelings of r„ for which we were able to compute the partition functions: we refer to 
them as the "Schreier" labeling and the "rotation-invariant" labeling. 
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4.1. The "Schreier" labeling. In the "Schreier" labeling, at a given corner of labeled r„ we have a copy 
of labeled r„_i reflected with respect to the bisector of the corresponding angle, see the picture below. 




It turns out that this labeling of the graph r„ can be alternatively described by considering the labeled 
Schreier graph I]„ of the Hanoi Towers group and then contracting the edges connecting copies of S]„_i 
in E„, hence the name "Schreier" labeling. 

For every n > 1, the number of vertices of r„ is |y(r„)| = |(3"~^ + 1). This implies that, for n odd, 
|y(r„)| is odd and so we allow dimer coverings touching cither two or none of the corners. If n is even, 
|y(r„)| is even, and so we allow dimer coverings touching either one or three corners. We say that a 
dimer covering of r„ is: 

• of type /, if it covers no corner; 

• of type g°'^ (respectively g"^ , g^"^) , if it covers the corner of r„ where two edges a and h (respectively 
a and c, h and c) meet, but does not cover any other corner; 

• of type h"'^ (respectively /i"^, h^'^), if it does not cover the corner of r„ where two edges a and h 
(respectively a and c, b and c) meet, but it covers the remaining two corners; 

• of type t, if it covers all three corners. 

Observe that for n odd we can only have configurations of type / and ft., and for n even we can only have 
configurations of type g and t. 

From now on, wc will denote by fn,gn^, g^f, g^, h'^, h^^f, ft^^, t„ the summand in the partition function 
<I>„(a, 6,c) counting the coverings of the corresponding type. For instance, for n = 1, the only nonzero 
terms are /i = 1, ft"'' = c, ft°'^ = b, = a, so that we get 

$i(a, 5, c) = 1 + a + 5 + c. 

For n = 2, the only nonzero terms are t2 = 2abc, g2^ = 2ab. g'^'^ — 2ac, g\'^ ~ 2bc, so that we get 

$2(0, b, c) = 2{abc + ab + ac + be). 

In the following pictures, the dark bullet next to a vertex means that this vertex is covered by a dimer in 
that configuration. Since the graph r„+i consists of three copies of r„, a dimer covering of r„-|_i can be 
constructed from three coverings of r„. For example, a configuration of type / for r„+i can be obtained 
using three configurations of r„ of type g""^ , g"^"^ , g'"^ , as the first of the following pictures shows. 
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A 

AA 




A 



A 

A/ 




/\A A/\ 



By using these recursions and arguments similar to those in the proof of Theorem l3.1l one can show that 
for n odd (respectively, for n even), the number of coverings of types /, /i"'', ft."'^, h^'^ (respectively of types 
i, 5"^ s'"') is the same, and so is equal to the quarter of the total number of dimer coverings of r„. 

Theorem 4.1. For each n, the partition function of the dimer model on r„ is: 

( 3" -1-3 

I $„(a, 6, c) = 2(4a6c) J {abc + ah + ac + he) for n even 



$„(a,6, c) = (4a6c)" 



■{l + a + h + c) 



for n odd. 



Hence, the number of dimer coverings of r„ is 2 2 



Proof. The recursion shows that, for n odd, one has: 



(6) 



Similarly, for n even, one has: 



fn - 




ac 
itin- 


bc 
li/ri-1 


uab 
"n 


— 2f 


-l9n- 


be 
itin- 


uac 
"n 


= 2tn- 


^ab 


bc 
-l9n- 


ubc 
^"-11 


— It 


^ab 

-l9n- 


ac 
ifJn- 



(7) 



^ab 



^bc 



nuab uo-c ubc 
^"'n-l"n-l"n-l 

- ^Jn-l"'ri-l"-n-l 

_ O f uab t,bc 

- ^Jn-l"-„-l"-„-l 

-Of h*^^ 

- ^^n-l"-Tt-l"-n-l 



The solutions of systems ^ and ([7]), with initial conditions /i 
computed by induction on n. We find: 



= 1,/if = c,/ij'= = 6,/i!;^ = a, can be 



tn = 2abc{Aabc) 
gl^ = 2ah{Aahc)- 
g^," = 2ac{4:abc)- 
= 2hc{4ahc)- 



n even. 



The assertion follows from the fact that $„(a, &, c) = /„ + /i„ + + /i„ for n odd and $„(a, 6, c) 



/„ = (Aahc) — 
hf = c(4a6c)- 
= 6(4a&c)- 
= a(4a6c)- 



n odd. 



^71 



9' 



ab 



'9n 



g„ for n even. The number of dimer coverings of r„ is obtained as $n(l, 1, !)■ 



□ 
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Corollary 4.2. The thermodynamic limit is ^log(4a6c). In particular, the entropy of absorption of 
diatomic molecules per site is ilog2. 

The number of dimer coverings and the value of the entropy have aheady appeared in [3] , where the 
dimers on r„ with the weight function constant 1 were considered. 

Note also that the number of dimer coverings found for Sierpihski graphs r„ coincides with the number 
of dimer coverings for the Schreier graphs I]„ of the group i?*-"^^ (see Section [S]). 



4.2. "Rotation-invariant" labeling. 




This labeling is, for n > 2, invariant under rotation by For n > 3, the copy of r„_i at the left 
(respectively upper, right) corner of r„ is rotated by (respectively 27r/3, An/ 3). 




We distinguish here the following types of dimers coverings: we say that a dimer covering is of type g 
(respectively of type h, of type /, of type t) if exactly one (respectively exactly two, none, all) of the 
three corners of r„ is (are) covered. Observe that by symmetry of the labeling we do not need to define 
gab^ gac^ gbc^ l^ab^ j^ac ^ j^bc _ gg^gy chcck that in this modcl for n odd we can only have configurations 
of type / and h, for n even we can only have configurations of type g and t. By using recursion, one also 
checks that for n even, the number of coverings of type t is one third of the number of coverings of type 
g, and that for n odd, the number of coverings of type / is one third of the number of dimer coverings of 
type h. 

Next, we compute the partition function $„(a, 6,c) associated with the "rotation-invariant" labeling. 
We will denote by fm gm hn, tn the summands in (a, 6, c) corresponding to coverings of types /, g, h, t. 
For instance, for n = 2, we have g2 = 3c(a + b), t2 ^ a'^ + 6"^, so that 

$2(a, b, c) = a^ + b^ + 3c(a + b). 

Theorem 4.3. The partition function ofT„, for each n>2, is given by: 

$„(a, &, c) = 2 2 (a3 -I- 3 j-gj, _|_ ^^y j ^^^3 ^3 _|_ 3^^^, _|_ ^^-^ j^^^ ^ ^^^^ 

$„(a,6,c) = 2'^ 2"'(a3 + &3)^ {ac + be)'' i'\3{a^ + b^) + c{a + b)) for n odd. 



Proof. Similarly to how we computed the partition function for the "Schreier" labeling, we get, for n > 3 
odd: 

(8) /-^""^^(^V") 
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and for n even: 



(9) 



9n = 6/„_i (^3^^ 



The solutions of systems ([5]) and with the initial conditions 92 = 3c(a + 6) and t2 = + b^, can be 
computed by induction on n: one gets, for n even, 

t„ = — (a^ + — 4 — (ac + &c) — J — 
.g„ = 3 • 2^^(a3 + 63)^^^(ac + 6c)^^ 

and for n > 3 odd 

f /„ = 2'-^^ (a3 + 63) (ac + 6c) 

\ /i„ = 3 • 2^^^ (a3 + 63) (ac + 6c) . 

The assertion follows from the fact that $„(a, 6, c) = /„ + /7,„ for n odd and $„(a, 6, c) = t„ + gn for n 
even. □ 

Corollary 4.4. T/ie thermodynamic limit is i log 2 + -j^ log(a3 + 63) + i log c(a + 6). 
By putting a = 6 = c = 1, one finds the same value of entropy as in the "Schreier" labeling, as expected. 

5. Statistics 

In this section we study the statistics of occurrences of edges with a given label in a random dimer 
covering, for the Schreier graphs of _ff^3) ^j^^^^ for the Sicrpihski triangles. We compute the mean and 
the variance and, in some cases, we are able to find the asymptotic behavior of the moment generating 
function of the associated normalized random variable. 

5.1. Schreier graphs of the Hanoi Towers group. Denote by c„ (by symmetry, a„ and 6„ can be 

studied in the same way) the random variable that counts the number of occurrences of edges labeled c 
in a random dimer covering on I]„. In order to study it, introduce the function $5i(c) — $^(1, l,c), for 
i = I, II, III, IV, and observe that $^^^(c) = $^^(c). Moreover, denote by ^n.i and cr^j the mean and 
the variance of c„ in a random dimer covering of type i, respectively. Note that we have finjii = f-n,iv 
and aljjj = aljy. 

Theorem 5.1. For each n> 1, 

3"-i + l 3"-i+3 3"-i-l 

MnJ — ^ , t^nJI = ^ , f-nJII = ^ , 

2 _ 3" - 6n + 3 2 _ 3" + lOn - 13 2 _ 3" - 2?i - 1 

Proof. For a = b = 1 , the system ^ reduces to 

^11 ^ KO! f^^")^ 
= /$///y3 , KK'K" 

with initial conditions ^{(c) = c, ^{^(c) = c^ and $(^^(c) = 1. Now put, for every n > 1, q„ = and 
fn — £1 ■ Observe that both qn and r„ are functions of the only variable c. In particular, for each n. 
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one has 9n(l) = ^ri(l) — 1, since the number of dimer covering is the same for each type of configuration. 
By computing the quotient <i>^^i/$^_|_i and dividing each term by ($^J'^, one gets 



(10) Qn+i = T^-^. 

7 + Inri 



Similarly, one has 

(11) rn+i 



' n ~ c 

7 + Inrl 



Using ((To)) and (fTT|) . one can show by induction that g^j(l) = 1 and ^^(l) = —1, for every n > 1. 
Moreover, (j"(l) = 4(n - 1) and r^(l) = n+l. 

From the first equation of the system, one gets $fj(c) — "~ (l + g„_i(c)r^_2 (c)) . By applying the 
logarithm and using recursion, we have: 

n— 2 n— 1 

log($^(c)) =3"-ilog($((c)) - ^3Mogc + ^3"-i-nog(l + cg,(c)r2(c)). 

k=0 k=l 

Taking the derivative and putting c = 1 , one gets 



f^nj — 



_ 3"-i + 1 

c=l — 7^ : 



what is one third of the total number of edges involved in such a covering, as it was to be expected 
because of the symmetry of the graph and of the labeling. Taking once more derivative, one gets 

C(c)$^(c)-«(c))^ _ _ 3"-^-6n + l 

Hence, 

2 _ <(1) , Ka)? ^ 3"-6n + 3 

"■^ $,^,(1) ^ mw mw 4 

In a similar way one can find jinji, '^nji^ l^^njii) ///■ 

Observe that one has /^„j/ > /i„j > fJ.n,iii- this corresponds to the fact that the distribution of labels 
o, b, c is uniform in a configuration of type I, but not in the other ones. In fact, a configuration of type 
II has a loop labeled c, but a configuration of type III (resp. IV) has a loop labeled b (resp. a): so the 
label c is "dominant" in type II, whereas the label b (resp. a) is "dominant" in type III (resp. IV). 

5.2. Sierpiiiski triangles. 

Theorem 5.2. For Sierpinski triangles with the "Schreier" labeling, for each n > 1, the random variable 
Cn has 

3"-^ 2 3 

- ^> - Ye- 

Moreover, the associated probability density function is 

3s{x + ^) + \S{x-V3) n odd 



where d denotes the Dirac function. 
Proof. Putting a = b = 1, one gets 



jS{x — + jS{x + -v/S) n even, 



$„(c) = (4c) J (c + 3) for n odd 
^$„(c) = 2(4c)— r^(3c+ 1) for n even. 
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The mean and the variance of c„ can be computed as in the previous theorem, by using logarithmic 
derivatives. Now let C„ = be the normalized random variable, then the moment generating 

function of C„ is given by 



We get 

( ,/^. „ 

n odd 



n even. 



and the claim follows. □ 



Observe that the moment generating functions that we have found only depend on the parity of n. The 
following theorem gives an interpretation of the probability density functions given in Theorem 15.21 

Theorem 5.3. For n odd, the normalized random variable Cn is equal to -y/S in each covering of type 
h'^^ and to ^ in each covering of type /, /i""^, /i^"^. For n even, the normalized random variable Cn is 

equal to — a/S in each covering of type g"^^ and to -y= in each covering of type t , g"''^ , g'"^ . 



Proof. The assertion can be proved by induction. For n = 1,2 a direct computation shows that the 
assertion is true. We give here only the proof for n > 2 odd. The following pictures show how to get a 
labeled dimer covering for r„, n odd, starting from three dimer coverings of r„_i. One can easily check 
that these recursions hold, by using the definition of the labeling of r„. 
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If we look at the first three pictures, we see that the variable C„ in a dimer covering of type /, /i"'^, h^'^ 
is given, by induction, by the sum of two contributions 1/^/3 and one contribution — -s/S, which gives 
— 1/^/3. The fourth picture shows that the variable C„ in a dimer covering of type /i"'' is given, by 
induction, by the sum of three contributions I/a/S, which gives -\/3- A similar proof can be given for n 
even. The statement follows. □ 

Theorem 5.4. For Sierpinski graphs with the "rotation-invariant" labeling, for each n > 2, the random 
variables a„ and 6„ have 

3"-i 2 4 • 3"-i + 3 
A^n = ^ <yn = ^ 

and they are asymptotically normal. The random variable c„ has 



and the associated probability density function is 



\5{^x — -^) + \5{^x + \/3) for n even 
l5{x + 4^) + \5{x - V3) for n odd. 



Proof. By putting 6 = c = 1 in the partition functions given in Theorem 14.31 one gets 



j $„(a) = 2'' 2"'(a3 + 1)^—4-^ (a + 1)'' [a^ + 3a + 4) for n even 
Similarly, one can find 



i-2_ 



+ 4) for n odd. 



1-1 3"-i-3 



I $„(c) — 2 2 c I (3c + 1) for n even 
[<i)„(c)=:2 2 c I (c + 3) for n odd. 
Then one proceeds as in the previously studied cases. □ 

A similar interpretation as in the case of the "Schreier" labeling can be given. 

Theorem 5.5. For n even, the normalized random variable Cn is equal to —VS in each covering of type 
t and to in each covering of type g. For n odd, the normalized random variable C„ is equal to VS in 

each covering of type f and to in each covering of type h. 

Remark 5.6. In [TB] the authors study the statistical properties of the dimer model on r„ endowed with 
the "directional" labeling: for n even (which is the only case allowing a perfect matching), they get the 
following expressions for the mean and the variance of the number of labels c: 

3"-i + l 2 3"-i-3 

Mn = ^ = 1 • 

Moreover, they show that the associated normalized random variable tends weakly to the normal distri- 
bution. 



COUNTING DIMER COVERINGS ON SELF-SIMILAR SCHREIER GRAPHS 



29 



References 

[1] L. Bartholdi and R. I. Grigorchuk, On the spectrum of Hccke type operators related to some fractal groups, Tr. Mat. 
Inst. Steklova 231 (2000), Din. Sist., Avtom. i Beskon. Gruppy, 5-45; translation in Proc. Steklov Inst. Math. (2000), 
no. 4 (231), 1-41. 

[2] L. Bartholdi, R. I. Grigorchuk and V. Nekrashevych, From fractal groups to fractal sets, in: "Fractals in Graz" (P. 

Grabner and W. Woess editors). Trends in Mathematics, Birkauser Verlag, Basel, 2003, 25—118. 
[3] I. Bondarenko, Groups generated by bounded automata and their Schreier graphs, PhD Thesis Texas A&M, 2007, 

available at http://txspace.tamu.edu/bltstreaiii/handle/1969. l/85845/Bondarenko.pdf?sequence=l 
[4] S.-C. Chang and L.-C. Chen, Dimer coverings on the Sicrpinski gasket, J. Stat. Phys. 131 (2008), no. 4, 631—650. 
[5] D. D'Angeli, A. Donno and T. Nagnibeda, Partition functions of the Ising model on some self-similar Schreier graphs, 

in: "Random Walks, Boundaries and Spectra" (D.Lenz, F. Sobieczky and W. Woess editors). Progress in Prob. 64 

Birkhiiuser, Springer Basel, 2011, 277-304. 
[6] D. D'Angeli, A. Donno, M. Matter and T. Nagnibeda, Schreier graphs of the Basilica group, J. Mod. Dyn. 4 (2010), 

no. 1, 167-205. 

[7] M. E. Fisher, On the dimer solution of planar Ising models, J. Math. Phys. 7 (1966), no. 10, 1776-1781. 
[8] R. I. Grigorchuk, Solved and unsolved problems around one group, in: "Infinite groups: geometric, combinatorial 
and dynamical aspects" (L. Bartholdi, T. Ceccherini-Silberstein, T. Smirnova-Nagnibeda and A. Zuk editors), Progr. 
Math., 248, Birkhiiuser, Basel, 2005, 117-218. 
[9] R. I. Grigorchuk and Z. Sunic, Self-similarity and branching in group theory, in: "Groups St. Andrews 2005, I" , London 
Math. Soc. Lecture Note Ser., 339, Cambridge Univ. Press, Cambridge, 2007, 36-95. 
[10] R. I. Grigorchuk and A. Zuk, On a torsion-free weakly branch group defined by a three-state automaton. International 
Conference on Geometric and Combinatorial Methods in Group Theory and Semigroup Theory (Lincoln, NE, 2000). 
Intemat. J. Algebra Comput. 12 (2002), no. 1-2, 223-246. 
[11] M. Gromov, Structures metriques pour les varietes riemanniennes. Edited by J. Lafontaine and P. Pansu, Textes 

Mathematiques, 1. CEDIC, Paris, 1981. 
[12] P. W. Kasteleyn, The statistics of dimers on a lattice I. The number of dimer arrangements on a quadratic lattice, 
Physica 27 (1961), 1209-1225. 

[13] P. W. Kasteleyn, Graph theory and crystal physics, in: "Graph Theory and Theoretical Physics" (F. Harary editor). 

Academic Press, London, 1967, 43-110. 
[14] R. Kenyon, Lectures on dimers. Statistical mechanics, 191-230, IAS/Park City Math. Ser., 

16, Amer. Math. Soc, Providence, RI, 2009. Lecture notes for a course, available at 

http : //www . math . brown . edu/~rkenyon/papers/dimerlecturenotes . pdf 
[15] R. Kenyon, A. Okounkov, S. Sheffield, Dimers and amoebas, Annals of Math. (2) 163 (2006), no. 3, 1019-1056. 
[16] V. Nekrashevych, Self-similar Groups, Mathematical Surveys and Monographs, 117. American Mathematical Society, 

Providence, RI, 2005. 

[17] H.N.V. Temperley, M. Fisher, Dimer problem in statistical mechanics - an exact result. Philosophical Magazine 6 
(1961), no. 68, 1061-1063. 

[18] E. Tcufl and S. Wagner, Exact and asymptotic enumeration of perfect matchings in self-similar graphs, Discrete Math. 

309 (2009), no. 23-24, 6612-6625. 
[19] E. Tcufl and S. Wagner, Enumeration of matchings in families of self-similar graphs. Discrete Appl. Math. 158 (2010), 

no. 14, 1524-1535. 

Departamento de Matematicas, Universidad de LOS Andes, Cra 1 n. 18A-12 Bogota, Colombia 
E-mail address: dangeliOuniandes . edu. co 

Dipartimento di Matematica, Sapienza Universita di Roma, Piazzale A. Moro, 5 00185 Roma, Italia 
E-mail address: donnoOmat.uniromal.it 

Section de Mathematiques, Universite de Geneve, 2-4, Rue du Lievre, Case Postale 64 1211 Geneve 4, Suisse 
E-mail address: Tatiana. Smirnova-NagnibedaOunige . ch 



